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In order to close the set of relativistic hydrodynamic equations, constitutive relations for the 
dissipative fluxes are required. In this work we outline the calculation of the corresponding closure 
for the heat flux in terms of the gradients of the independent scalar state variables: number density 
and temperature. The results are compared with the ones obtained using Marie's approximation in 
which a relativistic correction factor is included for the relaxation time parameter. It is shown how 
including such correction this BGK-like model yields good results in the relativistic case and not 
only in the non-relativistic limit as was previously found by other authors. 
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I. INTRODUCTION 



The interest in relativistic non-equilibrium thermodynamics has grown recently, mostly due to recent experiments as 
well as new ideas in theoretical modeling [H-Q • Most works dealing with applications of the theory do not refer to the 
corresponding system of transport equations including dissipation to first order in the gradients. This is principally 
due to the fact that for decades such system has been thought to predict unphysical dynamics in the fluid However, 
it was recently pointed out that when kinetic theory is used to formulate the constitutive equations that complete 
such set, instead of phcnomcnological arguments, the system of transport equations is free of instabilities and shows 
no causality issues in the linear regime Q. 

Relativistic kinetic theory to first order in the gradients has been vastly explored by previous authors • In those 
works the purely relativistic term that appears in the closure for the heat flux is included in a generalized thermal 
force with one transport coefficient i.e. a relativistic thermal conductivity. However, the system of transport equations 
constitutes a set of five balance equations and needs to be eventually expressed in terms of only five state variables. 
In the usual representation, these variables are the number density n, the hydrodynamic four velocity W , which only 
has three independent components, and temperature T. Because of this, constitutive equations for the dissipative 
fluxes in terms of gradients of those quantities are required to close the system. The aim of this work is to show 
the highlights of the calculation of the heat flux from the complete Boltzmann equation using the Chapman-Enskog 
approximation within the n, IA V , T representation. The details and further analysis of this calculation can be found 
in Ref. fioj . In order to perform this task we consider the density and temperature gradients as independent forces 
and write two uncoupled integral equations. The coefficients involved in the constitutive relation are calculated and 
plotted for a constant cross section model. Here, we also compare the results with those found in previous works using 
a BGK-like approximation 0, [uj. In order to make such comparison, we adjust the parameter included in Marie's 
kernel by a relativistic factor arising from the structure of the kinetic equation. We further show that when this factor 
is taken into account, the deviation from the exact values are as significative as in the non-relativistic case. 

The rest of this work is organized as follows. In the second section some basic elements of relativistic kinetic theory 
are briefly reviewed including the Chapman-Enskog expansion for the solution of the Boltzmann equation. In the 
third section an outline of the calculation of the heat flux using two independent forces is shown. The results for the 
transport coefficients are compared to the values obtained in the BGK approximation in section four and some final 
remarks and conclusions arc included in the last section. 



II. RELATIVISTIC KINETIC THEORY 



In the kinetic theory of gases, the evolution of the distribution function of a a dilute gas is given by the Boltzmann 
equation. If the temperature T of such gas is high enough such that the relativistic parameter z = fcT/mc 2 , where m 
is the rest mass of the particles, c the speed of light and k Boltzmann's constant, is not negligible, the non-relativistic 
description is not appropriate. Instead, relativistic kinetic theory (RKT) constitutes a more accurate framework. In 
RKT the fact that the speed of the molecules can be close to the speed of light is taken into account. The evolution 
of the distribution function is then given by the special relativistic Boltzmann equation [l2| 

v a .f, a = j(.f.n, (i) 

where v a = 7 (a/, c) is the molecular four-velocity, 7 being Lorentz factor , and J(f, /') is the collision kernel @ 



J (//') = / / [f'fi'-ffi]Fo-(si)cindvi 



In the previous equation T the invariant flux and a d£l is the differential cross section H . The system here 
considered is a dilute, single component, neutral, non-degenerate, inert gas. The local variables, number density n, 
hydrodynamic four-velocity u u and internal energy density e, are defined as 



fldv*, (2) 
nu v = J fv u dv*, (3) 
ne = mc 2 / f~/ 2 dv*. (4) 
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respectively. The local equilibrium distribution function, which is the solution of the homogeneous Boltzmann equa- 
tion, in the relativistic case is given by Ref. [HI 



(o) 



1 



4ttc 3 zK 2 (I) 



exp 



(5) 



where K n (x) is the modified Bcssel function of the second kind. As shown in Ref. [Tj] , the heat flux can be calculated 
as the average of the chaotic kinetic energy, that is 



cf = mc 2 h£ / jv u fd*v 



(6) 

The chaotic velocity corresponds to the velocity measured in a local comoving frame. Thus, one can replace the 
invariant u v v v in Eq. ([5]) by its value in the local comoving frame, that is u v v v = —c 2 ~f. In this context the total 
energy and momentum fluxes for an arbitrary observer are given by the Lorentz transformed dissipative fluxes as 
calculated in the local comoving frame. This argument is clearly discussed in Refs. [Tolflil]. 

In order to calculate the heat flux, given in Eq. (|6|), one does not need an explicit solution to the Boltzmann 
equation. Using Hilbert's method one can obtain a good approximation to the integral in Eq. ©. Following the 
prescription of such method, the velocity distribution function is written as / = /'°) + /W where /W is considered 
to be a first order correction in the gradients of the state variables. Also, in order to associate the local variables to 
the equilibrium state, the following subsidiary conditions need to be imposed 



(i) 



dv* 



0. 



Under these hypotheses, the linearized relativistic Boltzmann equation, obtained by substituting / in Eq. {TJ 
considering only linear deviations from local equilibrium, reads 

WO) = /( o) c (0) 



(7) 



and 



(8) 



where 



C {<t>) = / / Wi + <t>' -<f>i- <£} A (0) ^ (fi) dVLdvl 



is the linearized collision kernel. Substituting on the left hand side and using the identity V a A >a = W v v v A n 



v v A^ Eq. (HJ yields 



vaf (0) = f{ 0) htl » v (n^ 



/(0) 



7 C { 7i - T 



lit. 

T 
1 



_1_ 

2: 



K 1 g)+K 3 (k) - 



(9) 



Notice that the first line in Eq. only includes the spatial components, since in the comoving frame h^ 4 = while 
the time derivatives are isolated in the second line which we have simplified using that 



-v u = -u"7<S£ 



This separation between temporal and spatial parts is particularly useful since for the Chapman-Enskog solution to 
first order in the gradients to exist, the previous order equations have to be satisfied. Because of this, the relativistic 
Euler equations namely, 



-nu. r 



ne p \ 

c 2 + c 2 I P^ n 0' 



(10) 



T 



Tf3 
%C n k 



a-, 
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are substituted in the expression in Eq. ((9j) such that only spatial gradients appear in it. Also, since the momentum 
balance equation is in terms of the gradient of the hydrostatic pressure, we introduce in it the relation 

dp\ (dp 



in order to be consistent with our representation. After Euler Eqs. (|10[) . arc introduced in Eq. Q and using the 
equation state for an ideal gas p = nkT, we can write 

Since in this work we are only interested in the heat flux, we have ignored the term proportional to ur a in Eq. (|12[) 
due to Curie's principle. 

Following Hilbert's method, the general solution for the integral equation (fT2"]) is given by the sum of a particular 
solution and a linear combination of the collision invariants as solutions to the homogeneous equation. That is 

T n 

= A(j)v a h a ^ + B{^)v a hf-^ + a + a v v v (13) 

where A("f) and #(7) are functions of 7 and the local variables n, v, T while a and a v are constants. By imposing the 
subsidiary conditions, required for the uniqueness of the solution, Eq. ([7]), 4> can be rewritten as 

cb = A( 1 )vPh^+B( 1 )v?h^. (14) 
1 ' n 

The coefficients .4(7) and B(-y) are further expressed in terms of orthogonal polynomials in the form 

00 

OO 

B(7) = £&*Cn(7). 

71=0 

with C n (j) = X)fc=o a knJ k - The elements of this set satisfy the orthogonality condition 

' £„( 7 )£m(7) exp"5 ( 7 2 - l) 3 / 2 rf 7 - S nm , 



More details on these polynomials can be consulted in Refs. [oL fTojl . Direct substitution of coefficients ^4(7) and 23(7) 
in the solution (fT4"]) leads us to the expression 

^ = h a J2 a n C n ( 7 y^ + h%J2 WtK — , (15) 

n=0 n=0 

It is worthwhile pointing out that the first-order approximation for the distribution function is obtained here in 
terms of the density and temperature gradients as was already discussed in d, [To[ • This is because of the chosen 
representation in which n and T are considered independent state variables. Thus, when the heat flux is calculated from 
the distribution function, the term proportional to the temperature gradient will lead us to a Fourier-like law, while 
the term proportional to the density gradient is a purely relativistic contribution. The latter has no non-relativistic 
counterpart since the corresponding term in Eq. (|12[) vanishes in such limit. 

Substituting the proposed solution 4> in the linearized Boltzmann equation (jT2"j) leads to an integral equation 



which can be separated into two independent equations as follows 
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v p h% 1-7 



Ml) 



h% V (7) A 



n=l 



As mentioned above, the explicit complete solution for /W is not required in order to calculate the heat flux. As 
will be shown in the next section only one term in the infinite sums appearing in the previous equations is enough in 
order to get a good approximation for . 



III. HEAT FLUX 



The heat flux is physically the average of a chaotic energy flux. Because of this, it is calculated in the local comoving 
frame, as is clearly stated in Eq. (j6|). Substituting f^ x > = f^<f> with 4> given by Eq. (fT5l) in Eq. ((6]) yields 



AkczK 2 {\) 



where the integrals in the bracket are given by 

r = 



A(j)v k — 



^v v exp * 



l<*y, 



(16) 



(17) 



T v - 

1 b — 



n 



7« exp" 



'* \J~i 2 - Id'). 



(18) 



In Ref. [l(| it is shown how, of all coefficients a n and b n in Eqs. (|17[) and (|T51) . only one of each is required in order to 
obtain a, good first approximation to 1% and Ig. After a cumbersome calculation, which can be found to some detail 
in Ref. [l[|, and assuming constant collisions cross-section (see one obtains 



-h" a 



Ln 



L, 



where the transport coefficients Lt and L n can be written as 

3 ckT 



Lt 



64 7T<7 



■Hz) 



3 ckT 

64 7T(T 



where 



Hz) 



* s (±); + + 2 ) *,(!)) 



(19) 



(20) 



(21) 



/a(*) 



A' 



A, 



5£ 



j 3 A, 



+ (^ + 2)A 2 (D) • 



(22) 



and C/(l/z) = K^(\ j z~) f K 2 {\ j z) . In order to compare these results with the ones obtained in the relaxation time 
approximation we will explore Marie's model in the next section and find appropriate expressions for the functions 
fi(z) and f 2 {z) in such context. 



IV. COMPARISON WITH THE BGK APPROXIMATION 



The BGK approximation consists in replacing the right hand side of the Boltzmann equation with a simplified 
kernel in which all the details of the collisions are substituted by a parameter. The kinetic equation within relaxation 
models reads 



v a f, a = ->C(f-f m ). 



(23) 
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where the parameter JC is assumed to be independent of . The transport coefficients obtained using this method will 
have a strong dependence on this parameter. Thus, in order to compare the results obtained in the previous section 
using the complete kernel with the ones obtained with this model we need an approximation for the parameter K.. 

To understand the role of this constant we follow the standard reasoning. We assume that the one particle distri- 
bution function does not depend on the spatial coordinates such that we can write 



dt 

from where one obtains 



7? = -£(/-/ (0) ) 



fit) = exp 



(-T)Kf/Wf)« 



(24) 



Equation (|24j) motivates the association of the quantity 7//C with a characteristic relaxation time r. In the non- 
relativistic case, the analog to Eq. (|24[) leads to assuming r ~ K. In Ref. the 7 factor is not included in JC which at 
the end leads to results for the thermal conductivity that deviate substantially from the complete solution value in the 
relativistic and ultra- relativistic scenarios. Only when 7 — > 1, in the non-relativistic limit, the results are consistent. 

Here we consider r ~ j/JC and propose a reasonable estimate for the parameter K. from known characteristic times 
related to the collisions in the system. In non-relativistic kinetic theory, the time between collisions r c for hard-sphere 
particles of cross section a is given by 

1 



Ann(j{g) 



where (g) denotes the mean value of the relative velocity of two particles. Usually, the relaxation time r appearing 
in BGK model, is considered proportional to r c and since (g) is of the order of the mean velocity (v) and/or the 
adiabatic sound speed v s = '/3m, we propose for this case 

7 1 



fC 4niraV(z) 

where V(z) is the characteristic speed proportional to (v) or v s . Finally to approximate JC we propose 

1 1 

JC AmraViz)^)' 



(25) 



(26) 



in analogy with the non-relativistic case. 

Using Marie's model, Sandoval et al. @ obtained expressions for the transport coefficients which can be written in 
same fashion as Eq. (|20|) where the functions fi(z) and f2{z) are replaced by 

, , 16c / K 2 \ f 1 / 1 \ 1 / 1 \ z 1 Anna 




jW(z) 





(28) 



where we have included the M's subscript to indicate Marie's model solution. 

The results of substituting expression (|2"5)) into (|27|l and (f2"5)) are shown in Fig. [U where we have considered 
V(s) = V2(v). As can be clearly seen from the plot, the deviation of the approximated solution from the exact 
one is not as significant as without the 7 factor, which is shown in Fig. [5J Moreover, the difference between exact 
and approximated solutions when the 7 factor is taken into account is not grater than to the one obtained in the 
non-relativistic case which is shown in Fig. [31 for the sake of comparison. 

With the information in hand, one can also infer the most reasonable expression for the characteristic velocity V{z) 
with the help of both models. To accomplish this task, we write 



1 L T 

JC Lmt 64iranc 



K-2 



5£ 



*»(!)) 



7 



I/2I relativistic 




Figure 1. Relativistic coefficients /1 and fi. The continuous line for the complete kernel and the dashed line for Marie's model 
including the 7 factor in the relaxation parameter. 



fi relativistic 





Figure 2. Relativistic coefficients /1 and fi. The continuous line for the complete kernel and the dashed line for Marie's model 
without taking into account the 7 factor in the relaxation parameter. 



and 



1 L r , 



£ + 50(*)-M*) S 



K L Mn Mncan z* (f K 3 (f ) + + 2) K 2 (f)) 



and consider the best approximation to V(z) to be the one that makes both coefficients Lt's or L n 's match perfectly. 
The results are shown in Fig. 4 together with the characteristic velocity given by Eq. (|2l)|) and v s . Note that V c (z) 
is the value one should give to the characteristic relative speed for the relaxation time model to match exactly the 
results obtained with the complete kernel. 



fi non— relativistic 




.5 2.0 



Figure 3. Non-relativistic case with the complete kernel (continuous line) and the BGK model (dashed line) . Only the thermal 
conductivity factor /1 is shown since fi corresponds to a purely relativistic correction. 
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0.0 0.5 1.0 1.5 2.0 



Figure 4. This plot shows a comparison of the different characteristic velocities V c (z) (long dashed line), {v)(z) (dashed line) 
and v s (z) (continuous line) . 



V. FINAL REMARKS 

The relativistic hydrodynamic equations at the Navier-Stokes level need closure relations in order to constitute a 
complete set. In particular, the constitutive equation for the heat flux can be written in terms of the gradients of 
the temperature and the density when such quantities are considered independent state variables. The corresponding 
coupling coefficients have been calculated before assuming Marie's model @, |Tl[ . Here we outlined the procedure 
for the exact calculation, within Chapman-Enskog's approximation, and compared the results in the case of constant 
cross-section collisions. 

We re-analyzed Marie's proposal as a relaxation model and found that the quantity related to a relaxation time 
is actually best approximated by {^)/K,. With this correction Marie's model proves to be a good approximation not 
only in very mildly relativistic systems. Additionally we compare the different values that can be assigned to the 
characteristic microscopic velocity when approximating the value for K. with the one inferred from the exact solution. 

We finally wish to point out the similarity in structure of the second term in Eq. (fH))) with the Dufour effect. 
This is only in the sense that such relation implies the existence of a heat flux even in isothcrmic conditions due to 
a gradient in density. However, the coupling between heat and the density gradient in relativistic gases cannot be 
properly referred to as a cross effect since there is no reciprocal (diffusive) flux. The Dufour effect is the heat flux due 
to a diffusive force that appears in more than one component systems and has a conjugate effect called Sorct effect. 
The role of the term here studied in the case of a binary mixture is still not clear and will be addressed in the future. 
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